Bosonic and fermionic transport phenomena of ultra-cold atoms in ID optical lattices 
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Using the micro-canonical picture of transport - a framework ideally suited to describe the dynamics of closed 
quantum systems such as ultra-cold atom experiments - we show that the exact dynamics of non-interacting 
fermions and bosons exhibit very different transport properties when the system is set out of equilibrium by 
removing the particles from half of the lattice. We find that fermions rapidly develop a finite quasi steady-state 
current reminiscent of electronic transport in nanoscale systems. This result is robust - it occurs with or without 
a harmonic confining potential and at zero or finite temperature. The zero-temperature bosonic current instead 
exhibits strong oscillatory behavior that decays into a steady-state of zero current only in the thermodynamic 
limit. These differences appear most strikingly in the different particle number fluctuations on half of the lattice 
as a consequence of the spin statistics. These predictions can be readily verified experimentally. 



PACS numbers: 72.10.-d, 67.10.Jn, 03.75.Mn 

Ultra-cold atoms loaded in various optical lattices help vi- 
sualize and gain a deeper understanding of many exciting 
quantum phenomena. For instance, Hanbury-Brown-Twiss 
types of experiments measuring spatial correlations of den- 
sity fluctuations in non-interacting bosons [ 1 1 and fermions 
111 clearly demonstrate bunching and anti-bunching effects, 
which are definite signatures of the underlying spin statis- 
tics. On the other hand, bosonic and fermionic Mott insulators 
13] m have been realized experimentally and the statistics do 
not seem to play an important role in their thermodynamic 
properties if magnetic phenomena due to multi-component 
fermions are not considered. It is then clear that the ques- 
tion of how and when statistics affect an atomic system is of 
great fundamental interest and worth further investigations. 

In addition, recent progress in experimental techniques of 
ultra-cold atoms allows us to probe not only static structures, 
but also the dynamics of atoms in optical lattices. There have 
been experiments ||5] |6l showing that interactions (via colli- 
sions of atoms) can affect the motion of an atomic cloud when 
released from a non-equilibrium position. It has also been 
shown that mass and entropy transport of bosons in optical 
lattices can be unexpectedly slow (7]. Adding to the excite- 
ment are the experiments of Ref [8J which show how spin 
dynamics and phase transitions can be simulated using ultra- 
cold atoms. The possibility of using these systems to sim- 
ulate electronic devices and developing atom-based devices 
has also opened another research direction named "atomtron- 
ics" |9|. The often parallel investigations of spin-statistics 
and transport have therefore gained momentum in the study 
of ultra-cold atoms in optical lattices. 

Here we explore the effect of the underlying statistics on 
atomic transport phenomena. Since the deviation from the 
equilibrium could be significant, the applicability of the Kubo 
formula may be limited. Moreover, the atomic cloud in a real 
experiment should be modeled practically as a closed sys- 
tem. This could invalidate the popular Landauer formalism 
ifTOl which connects the system to reservoirs to study its static 



transport properties. Instead, we base our theoretical studies 
on the micro-canonical formalism (MCF) of transport intro- 
duced in Ref. ifTTI to study particle dynamics in nanoscale 
systems (see also Refs. | fT2l [T3l ). Specifically, this formal- 
ism explores transport properties of a closed quantum system 
where the dynamics of a finite number of particles - whether 
interacting or not - are followed in time 1 13 1. MCF is thus par- 
ticularly suited for simulating transport in ultra-cold atomic 
systems as done in actual experiments. This formalism allows 
us to make specific predictions on the different roles fermion 
and boson statistics play during the dynamical evolution of 
atomic gases. 

Here we consider the following experimental situation that 
may be realized in ID optical lattices - even though general- 
izations to higher dimensions are possible We consider 
non-interacting particles whose initial state corresponds to a 
given filling n. For a one-band model with moderate lattice 
depth, if the atoms are non-interacting fermions (bosons) of 
a single species, the ground state is a band insulator (quasi- 
Bose-Einstein condensate with all atoms in the lowest energy 
state). Then we remove part of the atoms - say, those on the 
right half of the lattice. This can be realized by a sudden exci- 
tation with a focused laser beam so that those atoms leave the 
trap while the lattice remains the same. This creates a vacuum 
region which then induces a current. Figure [T] demonstrates 
this experimental setup. 

In the MCF, the system is initially prepared in the ground 
state of one Hamiltonian Hq and then it evolves according to 
a second Hamiltonian of which it is not an eigenstate 1 1 1 1 . 
We formulate the evolution in terms of the correlation matrix 
Cij{t) = (G5'o|cJ(t)cj(t)|GS'o), where c| (ci)denotes the cre- 
ation (annihilation) operator of a particle at site i and \GSq) 
denotes the ground state of Hq (compare to Ref. [12J where 
the single particle states were used). For a system of N sites, 
the current through the middle is / = -~d{NL{t)) /dt. Here 
iVi(t) = J^iU'i (i)ci(t) counts the number of particles on 
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Figure 1. Schematic plot of the set up for fermions that we simulate. 
The gray dots emphasize that atoms may be in a superposition of 
locations. 



the left half of the lattice. To measure the current, one may fol- 
low Ref. Q and prepare several identical systems, then takes 
images at different times. The current is the rate at which 
atoms are transferred from the left to the right. 

To determine the initial state, we consider a tight-binding 
hopping Hamiltonian of non-interacting particles, Hq = 
E(.j) cU, + Z.immuj^d^ij - N/2 - l/2f c]c,, on 
a finite, non-periodic lattice. Here, t is the hopping coeffi- 
cient and it determines the unit of time, to = h/t and m, u, 
d are the mass of atoms, the trap frequency of a background 
harmonic potential, and the lattice constant, respectively. We 
define Y = mu'^cP /i \N\\h Y - 0.001 as in Refs. f5 , 6|. For 
fermions of one species, s-wave scattering and the Pauli exclu- 
sion principle naturally forbid observable interactions among 
atoms. For bosons we assume the background scattering is 
negligible. We remark that no dissipation mechanism is con- 
sidered in our study, which is consistent with Refs. ||5]|6l- 

The system is filled up to the given filling factor n and the 
ground state is best described in the energy basis with the cre- 
ation (annihilation) operator b\ (bk)- Here the index k labels 
different energy eigenvalues and should not be confused with 
the momentum. The unitary transformation Cj = ^k{U)jkbk 
diagonalizes the Hamiltonian as Hh — X^aLi ^kb\bk- When 
y = 0, one has tk = -2t cos[fc7r/(iV + 1)] and U-jk = 

, TV. For i < the corre- 
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lation function at zero temperature (T = 0) is {b\bk') — Skk' 
for fermions and {b\bk') — NSkk'^ki for bosons, where 
k = 1 denotes the lowest energy state. For fermions of 
one species, at n = 1 the lowest band is completely filled, 
thus there are no spatial correlations for fermions due to 
Pauli exclusion principle, but there are ample spatial conge- 
lations for bosons. At finite T, {b\bk') = rikSkk', where 
Uk = (exp[(efc—/i)/A:sT]±l)^^ is the corresponding thermal 
distribution and ^ is the chemical potential. We set fc^ = 1. 

The correlation matrix in real space can be writ- 
ten down using the unitary transformation as Cij = 

J2k,k'i^'')ki{blbk')Ujk'- In Fig. |2|c) we show the density 
profile of the ground state with n = 0.5 initially and Y ~ 
0.001 at T = 0, which resembles the wedding-cake structure 



0.6 ^ 



(f ,n=1) 



(a) 



0.3 t 


-0.3 



(f,n=0.75) 



f,n=0.5) 




(f.n=0.25) 



IIMPHi|MpMf>hii»^ 

(b,n=1) i^^si2 

, I , I J 




0.9 
0.6 
^ 0.3 

-0.3 
-0.6 



% 4000 / 



512 



r\Aj'v^v- 



' T=0 Y=0.00i (d) 



\ T/t=4 



T/t=8 My=o.oi"^^ 

V T=0 







25 t/t 50 



75 



Figure 2. (Color online) (a) Currents in a uniform lattice for fermions 
(top) and bosons (bottom) with initial n = 1 at T = after the atoms 
on the right are removed. The fermionic currents for initial fillings 
n = 0.75,0.5,0.25 are also shown (labelled next to each curve). 
Here A'^ = 512. (b) The bosonic current for TV = 8192 with initial 
filling n — 1. (c) The density profile for fermions with A'^ = 512, 
n = 0.5, and Y — 0.001 at T = 0. (d) Fermionic currents in a lattice 
with a harmonic potential for n = 0.5 and A'^ = 512. The solid lines 
correspond to T/i = 0, 4, 8 with Y = 0.001 while the dashed line 
corresponds to T/t = with F = 0.01. Here Y = muj'^S /i. 



of bosonic Mott insulators |7|. At t = all atoms on the right 
half are blown off by a laser beam whose frequency is differ- 
ent from that of the laser generating the optical lattices so that 
the lattice remains intact. The consequence is that [t — 0) 
becomes zero if i or j is on the right half, but the submatrix of 
the left half (both i, j < N/2) is not affected. 
The evolution of is determined by 
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{U^)k^{U),k,Dkk'{0)e^"'-''''^'■ (1) 



Here I'fefe'(O) = J2i',j=iiU'')k'j'UikCij{t = 0) is the initial 
coiTelation matrix in energy basis after the atoms on the right 
half are removed. The current from the left to the right is then 



2ttm(cT 2(t)cAr/2+i(i)). 



(2) 



Although bosons and fermions obey different statistics, the 
expression for the respective current is the same. The behavior 
of the cuiTent, however, is qualitatively different. 

Figure |2|a) shows the cuiTents of fermions and bosons on 
a uniform lattice of = 512 with initial n = 1 at T = 
after the atoms on the right half are removed. The fermionic 
current reaches a plateau after a short initial transient time. 
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This constant current has a duration proportional to the system 
size, indicating that fermions will develop a true steady state 
of finite current in the thermodynamic limit. For the finite- 
system case, the current varies abruptly when the electrons 
traveling to the right reach the far boundary and come back, 
as is seen in Fig.|2ja). Here we will focus on the physics prior 
to this finite-size revival behavior The almost constant current 
found here is similar to the quasi steady-state current (QSSC) 
shown in Ref. [12!| where a step-function bias is applied to a 
half-filled lattice. As in that case, it shows that a quasi steady 
state can originate dynamically in the absense of interactions, 
solely due to the wave properties of the system ifTTI . More- 
over, QSSCs can develop in higher dimensional geometries, 
as shown in Ref. ||12|. 

Importantly, the fermionic QSSC is observable in the pres- 
ence of a reasonable harmonic trap potential for both T = 
and finite T |14| as shown in Fig. |2|d). We note that the 
length of the initial transient time during which the current 
oscillates significantly depends on the initial state configura- 
tion and how the system is set out of equilibrium. 

In contrast, the T — Q bosonic current exhibits rapid os- 
cillations with a period ss to (Al/Iavg ^ 1, where AI is 
the amplitude of the oscillations and lavg is the current av- 
eraged over a period of time). Therefore, a QSSC is not 
found fl 51 . However, as it is also evident from Fig. 
and (b) that the oscillations of the bosonic current decrease 
with time. For bosons with initial filling n = 1 on a uni- 
form lattice of size N, the correlation matrix in energy ba- 
sis after the atoms on the right half lattice are removed is 
Dkk'it = 0) = Nj2^^^^j^^^UaUj,U^kUjk' at T = 0. 
In the limit iV — >■ oo it approaches NF{p)F{p'), where 
F{p) = pcos(p7r/2)/[(l -p2)7r]. F{p = odd) = 0, except 
F{1) — 1/4: and F{p — even) is finite. 

The current from Eq. (|2| then becomes 



/ = 2i ^ sin 



2tt 



cos( 



pn 
N + 1 



) — cos( 



p'tt 
N + 1 



Un/2..pUn/2+i,p'NF{p)F{p') 



(3) 



In the final expression only terms with {p — = even} 
and {p' — l,p — even} contribute and their contributions are 
equal. We focus on the long-time limit and set t = {N/2)to. 
One can make other choices as long as t/tQ ^ 0{N/2) and 
obtain the same conclusions. We consider {p = l,p' = even} 
and the final result is simply twice the current of this part. Due 
to the rapid oscillatory part of the first term for large t/to, only 
< p' /{N + 1) < l/\fN may contribute finitely. The current 
then approaches 
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dpsin[7V(l - cos(p7r)] sin \ cx iV"^/^ 



Therefore in the thermodynamic limit, the long-time current 
decays to zero for bosons. This analysis is supported by our 
simulations shown in Fig.[2];a) and (b). 

Another interesting question is how the filling factor affects 
the QSSC for fermions. If one performs the experiments sep- 



arately for different initial filling factors n = 1, 0.75, 0.5, and 
0.25 on a uniform lattice using fermions, the QSSC does not 
exhibit linear dependence on n as shown in Fig. |2|a). This 
non-linear dependence of the QSSC versus the filling factor 
originates from the band structure, which results in differing 
amounts of correlations between sites on the lattice. As the 
filling gets lower, the initially remaining fermions on the left 
half lattice are more correlated because the constraint from 
Pauli exclusion principle is less restrictive. This correlation 
effect competes with the counting of contributions from dif- 
ferent energy states. However, the initial filling does not affect 
the qualitative behavior of the QSSC of fermions if it is not too 
different from n — 1. 

Before discussing other differences of the two dynamics, 
let us first check if they can be observed experimentally. For 
optical lattices with moderate depth, with Vb and Ei^ the lat- 
tice depth and recoil energy, respectively, a typical magnitude 
of the hopping parameter t can be estimated from the param- 
eters given in Ref. |16|. For Vo/i?_R = 10, t w 3 nK which 
corresponds to to = 2.5 ms. Therefore the oscillations of the 
bosonic current are not too difficult to detect. Moreover, the 
typical holding time for recent experiments can be of the order 
of one second |5 1 so the QSSC of fermions, which is of the or- 
der of about 50io ^ 0.13s for a lattice of = 512 sites with a 
reasonable harmonic potential, should also be observable for 
T < 8t according to Fig.^d). 

We now show that particle number fluctuations can distin- 
guish the underlying spin statistics when the system is dy- 
namically evolving. Let us evaluate the instantaneous fluc- 
tuations in the total number of particles on the left half lattice, 
AN"^. Experimentally AA^^ may be measured by repeating 
the preparation and measurement at a given time so the fluc- 
tuations could be inferred. By definition. 



AA^^ = {ND - {Nl)' 



Here (Nl) = EfJii^l) + ^E^^^ (nS,), where h 
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(4) 

4 c, 



and rii = {fii). This expression involves expectation values 
of four operators and one can decompose them using Wick's 
theorem. The decomposition, however, is sensitive to the un- 
derlying spin statistics. After some algebra, one has 

N/2 

^n-i(l - r. 

i=l i<j 



^"^fermions 
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^"^bosons 
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) + 2E|Q,f. (6) 

i<j 

In contrast to the current, the expressions for the fluctuations 
are different for fermions and bosons. In the following we 
focus on T = 0. 

For fermions, < ni, \cij \ < 1 due to Pauli exclusion prin- 
ciple. Clearly, the correlation terms |cijp reduce the number 
fluctuations. For the initial filling n — 1, right after the par- 
ticles on the right half lattice are removed the instantaneous 
number fluctuations on the left half are zero because Pauli ex- 
clusion principle forces one particle per site. In contrast, the 
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Figure 3. Particle number fluctuations on the left half lattice for (a) 
fermions and (b) bosons in a uniform lattice. Here n = 1 before 
the particles on the right half are removed and A*" — 512. The inset 
shows the bosonic fluctuations of a hypothetical case where the initial 
configuration is identical to that of the fermions. 



initial boson distribution does not have this constraint. Thus 
by removing particles on the right half lattice, one removes 
half of the total number of particles on average so initial num- 
ber fluctuations are present. Moreover, Icy p contributes pos- 
itively to the fluctuations of bosons. 

To better visualize these results we plot in Fig.|3]the quan- 
tity AiV^(i) — AN^{t = 0) for fermions and bosons in a 
uniform lattice. For fermions the fluctuations increase slowly 
after an initial transient time. This slow change in fluctua- 
tions corresponds to the QSSC in Fig.|2](a). The behavior for 
bosons is very different. The contribution from spatial corre- 
lation, |cij p, decreases with time so the fluctuations decrease 
accordingly. The initial state of bosons is a quasi condensate 
so there is strong coherence among particles on different sites. 
As the system evolves and more particles move to the right, 
the coherence decreases rapidly so the number fluctuations 
are reduced when compared to the initial number fluctuations. 
Since there are about N"^ terms in the contribution from spa- 
tial correlations, the observed change in AA^^ for bosons is 
several orders of magnitude larger than that for the fermions. 

There are important differences when comparing our results 
with the Hanbury-Brown-Twiss type of experiments [T] |2l, 
where spatial correlations of density fluctuations are mea- 
sured. While those experiments [ 1 2] utilize the underlying 
spin statistics for an equilibrium system but implement a dy- 
namical measurement (time-of-flight imaging), our discussion 
depends on the statistics of the initial distribution as well as 
the corresponding quantum dynamics. Thus the difference in 
AA^^ demonstrated here is a dynamic property. 

To further stress this point, we consider a hypothetical case 
for bosons where the initial configuration is assumed to be 
exactly the same as that of fermions, i.e., the lowest N energy 
states for n = 1 are occupied by one boson each before the 
particles on the right half are removed. Due to this selected 
initial correlation matrix, the bosonic current will be exactly 
the same as the fermionic current shown in Fig.[2|a) and will 
show a QSSC. However, the bosonic number fluctuations, as 
shown in the inset of Fig. [3|b), are very different from the 
fermionic ones shown on Fig. [3|a) so the underlying statistics 



can still be resolved. 

We conclude that non-equilibrium dynamics of non- 
interacting atomic gases already demonstrates unexpected 
phenomena. Ultra-cold atoms in optical lattices have the great 
potential of demonstrating a QSSC and its dependence on 
the spin statistics in non-interacting clean systems, which is 
difficult to realize for electronic systems due to Coulomb in- 
teractions. Since interactions between atoms are experimen- 
tally tunable, we anticipate that transport phenomena in the 
presence of interactions will also reveal interesting physics 
ITZl . The setup discussed here could also be of interest to the 
community of quench dynamics ifTSl . Therefore the micro- 
canonical picture of transport presented here could contribute 
significantly to our understanding of many fundamental as- 
pects of transport phenomena in ultra-cold atoms. 
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